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ON THE PLANE POTENTIAL PLOW PAST A 
LATTICE OF ARBITRARY AIRFOILS 
By i. E. Garrick 



SUMMARY 



The tv/o -dimensional, incoiripressible potential flow 
past a symmetrical lattice of airfoils of arbitrary shape 
is investigated, theoretically. The problem is treated 
by usual methods of conformal mapping in several- stages, 
one stage corresponding to the mapping of the framev/ork 
of the arbitrary line lattice and another significant 
stage corresponding to .the Theodorsen method for the 
mapping of the arbitrary single, wing profile into a circle, 
A particular feature in the theoretical treatment is the 
special handling of' the regions at an infinite distance 
in front, of and behind the lattice. Expressions are 
given for -evaluation of the velocity and pressure distri- 
bution at the airfoil boundary. An illustrative numerical 
example is included . 



INTRODUCTION 



This paper treats the problem of determining the flow 
pattern, or the velocity and pressure fields, associated 
with the uniform flow past an infinite rov/ of symmetrically 
placed airfoils of the same shape. This airfoil-lattice 
problem occurs in the design of turbine blades, wind-tunnel 
vanes or grids, and elsev/here . There is a purely mathe^ 
matlcal interest in the problem that concerns the field 
of conformal mapping of infinitely connected regions. 
Analogous two-dimensional lattice" problems occur in the 
steady flow of heat and electricity. 

Considerable Ingenuity has been devoted to the 
airfoil-lattice problem, especially in the turbomachine 
studies in the German literature and more recently in the 
British studies; nevertheless, a survey of the available 
literature Indicates that nearly all the treatments 
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employed and the results obtained are of a special or 
indirect nature v/hich involve, for example, lattices of 
thin lines or approximate graphical procedures. Re- 
cently, however, A. Hov/ell in a British paper of 
limited circulation has v/ritten briefly on the theory of 
arbitrary airfoils in cascade, fiowell applies a special 
transformation to an airfoil lattice to convert the lattice 
region to a somewhat random, simply connected region and, 
with the. aid of several stages of conformal mapping, ob- 
tains a region about a circle. 

The problem of determining the incompressible poten- 
tial flow past an arbitrary single wing section v/as 
studied by Theodorsen (reference 1), who gave a practical 
procedure for its solution. The case of two wing sections, 
or the arbitrary biplane, was later treated in reference 2. 
The determination of the flow past an infinite lattice of 
airfoils of the same shape is a problem intermediate in 
difficulty in comparison v/ith the afore-mentioned ones. 
The treatment for resolving this problem given in the 
present report is similar to that for the arbitrary single 
wing section but the calculations are- more involved. 

The problem will herein be studied by the usual 
method of conformal mapping. It Is convenient to ac- 
complish the result in three or four stages: The airfoil 
lattice is first replaced by its skeleton, or framework 
. . .of line segments. The initial mapping function employed 
transforms the lattice skeleton into a circle. In the 
plane of this circle there are two singular points, known 
as .branch points. These points have dual significance: 
They correspond to infinite regions in front of and behind 
the lattice of lines, and they enter in the problem of 
reducing the lattice region (multiply connected region) 
to the region of a single body (simply connected region). 
If now an arbitrary airfoil shape is generated or given 
around the framev/ork of lines, then in the plane of the 
circle a circular-like contour is generated around the 
original circle. This contour may be transformed into 
an exact circle by the well-known procedure given in ref- 
erence 1 or 5. ■ The original tvifo significant -points are 
then traced by a transformation due to H. A, Schwarz. 
A final elementary transformation will bring the circle 
into a standard circle for which the two characteristic 
branch points are symmetrically placed. The region of 
this circle is considered the standard region for deter- 
mining the flovr pattern. 
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to 



For illustrative purposes an outline of a 'procedure 
for calculating pressure distributions is included. The 
method may be followed without reference to the theory 
by readers interested mainly in numerical results. For 
convenience, a list of syiiibols is given in appendix A. 



ANALYSIS 



Initial transformation for lattice of straight 
lines Consider the transformation (reference [|) 



log + log ^ 

b - z J 




(1) 



where g, b, and a are real nurabers and b > a. 
Introduce coordinates vj and 6 by means of the relation 



^i/+i9 

z ^ = ae 

and let 

b/a = e^o 

Equation (1) may then be expressed as 

cosh Yo cosh {\!/ + 19) 



cosh Yo - °osh + 18) 



(2) 



(5) 



If ij/ = 0, according to equation (2), z» lies on a 
circle^ of radius a (fig. 1(a)). According to equa- 
tion (4), i>-]_ = + iy^ is the logarithm of a real posi- 
tive function and consequently represents a real function 
(its principal value) and the infinite sequence of values 

differing from this function by •^:2kTTi,;, where k is 

any integer. The transformation illustrated in fig- 
ure Kb) is that of an infinite lattice of unstaggered 
lines of gap g in the t]^-plane into the circle of 



k 



radius a In the z^-^plane. The points z» = b and • 
z » = -b correspond to infinity in front of and behind ^ 

the lattice, respectively. The inverse points - tt 

2 . 
and z > = are inside the circle of radius a. 

b 

In order to introduce stagger, it is convenient to 
consider the transformation 



ih L b + z' . 
log 1 



2 \ b - 2 




v/here h is real. This transformation can be written 
with the use of equations (2) and (5) as 



sinh Yo - s^^^ 



(5) 



If \|/ = 0, the expression v/ithin the brackets is a com- 
plex number of unit magxaitude; hence, the logarithm is a 
pure Imaginary number plus an infinite sequence of numbers 
differing by ''Zftl. Then ^2 ^ ^2 ^^2 r^epresents a 

sequence of real numbers differing by h and the* lattice 
is one of horizontal lines displaced from each other by. 
h (fig. 1(c)). 

The transformation for the general staggered-llne 
lattice is a combination of equations (i|) and (5) 

t = ti + ^2 ■ (6a) 

or 



i = ± e-iP log ^-1^ log ^ I (6b) 

2tt \ b - z.'. 2, . af. 



a2 

z» + — 
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where 

gap g = d cos p 
stagger h = d sin p 
stagger ratio h/g = tan p 

the parameter d may be called the slant gap (fig. 1(d)), 
and p the stagger angle. 

The geometry of the lattice may be expressed in terms 
of the parameters Yq P noting that the chord 

length m.ay be obtained from the ( singular or critical) 
values of 6 which correspond to the end points of the 
chord and are solutions of the equation d^/dz' = 0. 
This equation gives the result 



tan 6 = tanh tan p 
' o ^ 



(7a) 



or, for later reference, 

g _ cosh Yp P 



cos 



sin e = 



sinh Yn P 



v/here 



Q = ^ccsh'^Y^ 



sm 



1/2 



Relations (7) ina^y be employed in two v/ays : (1) When the 
parameters Yq P ^'^^ given, the relation determines 

the two critical .values of 9, 6^. and 9^, where the 

subscripts I and t refer to leading edge and trailing 
pdge, respectively, and ®t ^ V/hen 9j or 

tan 9 J and the stagger angle p are given, the relation 

determines the parameter. Yq» 

The chord c may be obtained by putting 9 9^ 
and 9 = 8^' in equation (6a) and taking the difference 
in abscissas Xj and x^. Prom equations (ij.) to (7), 
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X 



t 



2d / = « -, Q + cos 3 ^ . - . -1 sin p\ 
= — I cos p log + sm p tan 1- ) 

" \ sinh Yo ^ / 



(8) 



By means of equation (8), the p^xpameter Yq ^® 
presented directly in terms of givon values of the gap- 
chord ratio for any staggo^r ratio, A representative 
chart relating gap-chord ratio, stagger angle, and 

is shown in figure 2; some values are given in table I. 

Inversion of equations ( I4.) to (6)*- The initial 
.transf onnations may be thotightof as raapping a framework 
of chords of an arbitrary lattice. into a circle. If a 
shape is generated around the chords in the 2* -plane, a 
contour is generated around the circle of radius a. 
This contour, which must exclude the points 2 ' ^ 

2' = b and must enclose the points z' = - — and 
2 ^ 
= may be considered to be completely defined by 

the function ^{Q) . If a lattice of airfoils is pre- 
assigned, the function \!/(6) must be found from the given 
coordinates of the airfoil shape. In order not to inter- 
rupt the sequence of main ideas, the details of this 
problem are relegated to appendix B, with certain remarks 
on the practical achievement of a nearly circular contour. 

Transformation of contour in -plane to circle in 
z -plane <>- It is assumed now that the circular-like contour 
in the -plane v/hich corresponds to the airfoil contour 
of the lattice is either given or determined; that is, 
the function ^{Q) is known in the boundary expression 
^J/+i0 

z ^ = ae . By the procedure of reference 1 or 5, the 
transformation 

z' = ze^^^^ (9a) 

where 

f(z) =f:£g = log^ (9b) 

1 z^^ z 

and Cri are complex coefficients determined by the 
boundary, is then employed to transf orra the 2 '-contour. 
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into a circle z = ae in the s-plane. The trans- 

forination (9a-) keeps the regions alike at infinity in the 
z^- and z-planes; that is^ z. = z* and dz*/dz = 1 at 
infinity. The corre.spondence of the boundaries is deter- 
mined by the functional equation 



cp -'9 = €((p) 





(10) 



for which a convenient numerical solution has been out- 

lined in reference 5- The radius of the circle R = ae'^ 
is determined by the relation 



M^) dcp (11) 



For consistency, the functional syrnbol ^(cp) is here used 
to denote the quantity ^ expressed as a function of cp - 

that is, ^r9(c?j( . In reference ^ the notations \!/{cp) 

and \l/j9tcp)j are used. 

It is necessary also to trace the correspondence of 
the points z^ =^ b. and z ^ =^ -b . Let z = p-|^ corre- 
spond to z^ = b and let z = correspond to z^ = -b. 
The values p-j^ and P2 ^^^7 determined by a relation 

(due to Schwarz) 'that expresses the value of a complex 
function in terms of an integral of the real part of the 
function along a circle. A simple derivation of the 
desired relation is shown in appendix C. The expression 
is 



log ^1 f(z) 
z 



^ xH<?) ^ (12) 

' z -iffi 
^ - R " 
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The values of ^-^ and inay be determined from . 

equation (12) by an iteration process that converges 
extremely rapidly. The. process may be described as 
follows: In equation (12), 1st the zeroth approximation 
to p-j^ be z = Zg = b and let the corresponding value 

of be written z» = Zq' = be"^^^^, v/here f(b) is 

the evaluation of equation (12) for z = b. It is 
actually desired, however, to have z' = b but, because 

Zt=:ZQf=:b + ZQ'-b 

the initial value of z' differs from the desired value 
by Zq^ - b. Furthermore, z = Zq differs from z = 

by approximately the same amount; hence, reducing zg 
by the quantity Zq ' - b gives 

Zi = Zq b - Zgt 

= b[2 . ef^(^)] 

which may be considered a first approximation to P^l* 

If it is desired to check this result or to obtain a 
second approximation, the process may be repeated; thus, 
from equation (12), find and 



Then, 



Zg = z^ + b - z-j^' 



which is a second approximation to p^^ and, in general, 
the nth approximation is 

^n ^ - ^n-l' 

It is clear that, should z^ correspond to ^n-^l^ ^n-l' 
must correspond to b' and the process automatically stops. 
The numerical process is given in appendix C; relatively 
elementary steps are Involved. In order to determine 
-p2, the process is applied v/ith b replaced by -b. 
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Transformation to standard circle In w-plane > - In 
order to obtain the flov/ pattern, it is desirable to 
introduce another . function v/hich traiisforras the circle 
in the z-plane into another circle in the v/-plane in. such 
a way that the characteristic points z = 3-[^' and 

z = map Into w = b and w = -rh, respectively. 

The region of the circle in the w-plane nay be considered 
the standard region. The desired transformation may be 
written as (see appendix D) 



(15) 




where 



K 



b^ + ^ 



ilk) 



and 

R = ae is the radius of the original circle in the 
z-plane, '^^ is the complex conjugate to ^j^, and S 

is the radius of the new circle in the v/-plane. The 
radius S is determined by 



S = be""'^l 



(15) 



where 



is obtained from 



cosh 



Yi = 



R 



r2 + 



^1 ^ 



(16) 



Complex velocity potential in v/-plane .~ Consider the 
flov/ function ^^(w) = $ + if , which is defined as 

o2 



V/ + 



ia T b + w ^ -ia ^ 
e log + e log 



w - ■ 



-Slog 



b - w 



(17) 
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The flov/ pattern may be regarded as due to a combination 
of singularities,, sinks, sources, and vortices, placed at 

the points w = tb and w = ±^ as indicated in figure 5 

It may be readily verified that the circle ot radius S - 

that is, w = Se"^^ - is part of a streamline and it may 
further be observed from figure 5 that the circulation 

around any contour which encloses the points w = 

and for which the points v/ = tb are exterior points, 

is r (positive if counterclockwise ) . The parameter a 

v/ill be interpreted "later as an angle of attack. 

The value of the circulation .V may be determined 

by means of the Kutta^ Joukowski condition for smooth flow 
at the trailing edge of the lattice. Let be the 

value of o on the boundary circle .Se-^^ that corre- 
sponds to the trailing edge of the lattice. The 
Kutta- Joukowskl condition then requires that the flow 
separate at 0 = o^, or that a stagnation point exist 
there . 

Vuth dQ/dw = 0 and w = Se'^^o, the following 
relation for P is found: 

r = - b sin (0^ + a) + ^ sin (a^ • a)] (l8) 

b2 - sit L ^ 

2 
b^ 

If s/b is replaced by e '"'l (equation (15)), equa- 
tion (l8) may be expressed as 

/cos sin \ 

T = -2Vd f ^ sin a + 2L cos a) (19) 

\cosh Y]^ sinh y-^ / 



Expr essions for ve l ocity in lattice field .- In order 
to obtain the flow pattern in the lattice field (t-plane), 
the component factors of the follov-.^ing expression are 
required: 

^ _ ^ dw dz dz ^ 
• dt ' dw dz dz ' dt 
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These terms may be obtained from equations (17)> (15)* 
(9), and (6). 

It is of particular interest to evaluate equation (20) 
explicitly for the regions at infinity in front of and 
behind the lattice and also on the lattice boundary itself. 
It is recalled that ^ =: co corresponds to = b, 

z = ^2_y w = b and that ^ = -co corresponds to = -b, 

z = -p2^- ^'^^ ~ combining terms according to equa- 

tion (20), the (reflected) inlet-velocity vector is 
obtained as 



dQ* 



= . Ve^^^^^^ - ^ eil 

2d 



(21) 



and the corresponding expression for the outlet-velocity 
vector is 



[da 



= - Ve ' e 

2d 



(22) 



By addition of equations (21) and (22), it becomes clear 
that the velocity vector of magnitude V and angle of 
attack a + 3 v/ith respect to the x-axis is one -half the 
vector sum of the inlet - and outlet velocities (fig. 1}.). 

If the angle of attack of the mean velocity vector 
v/ith respect to the- x-axis (chord direction) is denoted 
by = a + P, the velocity components in equations (21) 

and (22) are. 



Vx-|^ = -V cos + 



2d 



sin P 



Vy^ = -V sin + 



2d 



cos |3 
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and 



Vy^ = -V cos ay — =— sin P 

^2 .2d 



Vvo = V sin - — cos P 

X 2d 



The conventional angle of attack a is measured 
with respect to the noriiial to the slant line of the 
lattice. The components nomal to and along the slant 
line of the lattice, sometimes - referred to as "axial^^ 
and^ ^Uvhirl^* components, respectively, are found by rotating 
all' vectors in the xy-plane by angle p (fig. i|) . These 
components are, for the inlet velocity, 

Vj.j = -V cos a 



Vt = V sin a + — 
^1 2d 



and, for the outlet velocity, 

Vtj2 = cos a = Vfj]_ 

Vt^ = V sin a - ^ 
^2 2d 



The squares of the magnitudes of the inlet and outlet 
velocities are 



V2^ 



V2 



1 + 2 



2Vd 



sm 



a . M 

\zvdj 




v/here r/2Vd may be obtained from equation (19) • Observe 
that the inlet and outlet speeds are equal, V^^ = 

when a 0^ for any value of T. The inlet and outlet 
angles of attack with respect to the normal to the lattice 
line are 



13 



a]^ = tar 



sin a + 

,-1 2Vd 

cos a 



0^2 = tan 



sin a - 

-1 2Vd 

cos a 



and the angle through v;hich the streari is turned is 

2 -i— cos a 
-1 2Vd 



tan--^ — (2$) 

1 - 



a-t - ap — tan Tz 



^2Vd 

The component factors In equation (20) are now to 
be evaluated at the lattice boundary and, as the boundary, 
itself is part of a strefxinline, only the magnitudes of 
the factors are of interest. 

Prom equations (17) and (19) and with w = Se^*^, 

|sinh Yi sin a(cos a - cos Oq) 



(clQ 2Vd 1 



Idv/I- ttS cosh 2y]l " cios 2a 

+ cosh Yi ^os a(sin CJ - sin ^o)\ (2l|.) 

where the parameter Yi -^s defined in eqiiation (15)- 

In order to obtain dw/dz, it is convenient first 
to express equation (15) explicitly in w as 

b(l + K)z - b^Kp. - Pp) 

w = ^ (25a) 

(1 - K)z + KP^ + S2 

A standard form for the transformation of a circular 
region |zl= R into |w|= S is 



Ik 



z - 5 

w = RSe- -3 ^ (25b) 

R"- - 5z 

Comparison of equations (25a) and (25b) makes it clear 
that the complex parameter 5 and the real parameter X 
may be obtained from the following relations: 

6 = -Ll L2 (26a) 

1 + K ■ 

or, as a check relation, 

6 = "'"^ - 1' (26b) 

and 

or, by equating angles on both sides, 

X = argd + K) - arg^^8^ + p^) (2?) 

Prom equation (25b), the explicit correspondence of 

a point on the circle w = So"^*^ to a point on the circle 
i<p 

z - Re can be obtained as follows: 



gia ^ gi(cp+X) R 



1 - ,^ e-^^'' 



1 - Aei^ 



(28) 



R 

Let the complex number 6 be expressed as !5|e^^ and let 

1 - I e-i^ me^^ (29) 

where 

= 1 - 2 M cos ((P - T) + 1%. 

H j^2 



I 
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and 



|i(cp) = tan 



-1 



|5| 



sin (CD - T) 



1 - 



151 



cos (Cp - T) 



Observe that the denominator in equation (28) is the con- 
jugate of equation (29) and is therefore equal to me""^'^. 
There results for the correspondence of ■ a and cp 



a = cp + \ + 2p. 



(50) 



In particular, if the ( trailing-edge ) value of cp that 
corresponds to 9^ as determined bjr equations (7) is 

v/ritten as cp 



4. , 



v/nere 



'o t 

€(cp) at 9 = Qfc ^^'^^ equation (10), then 

^o =. ^o ^ + 2p.^ 
By differentiation of equation (25b), 
dw _ ns(R2 - 55) e^^ 



is the value of 



az 



(r2 - 

icp 



(51) 



On the boundary, put z = Re : then, the magnitude of 
equation is 



dw 


= -( 


dz 


R \ 



1 - 



p / 2 



(52) 



The expression for 



dz ' 



dz 



on the boundary is obtained 



from equation (9) in terms of the functions e(cp) and 
^(cp) of equation (10) as follows (see reference 3): 



dz' 
dz 



z 



1 + z 



df 
dz 



(55a) 



and, because f(a) on the boundary is 



f(z) = ^(cp) - % + i(e - cp) 
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where 



e - cp = c(cp) 



then 



! dz I z I \ dcp/ Vdcpy 



1/2 



(33b) 



The last factor of equation (20) is expressed from equa- 
tlon (6) on the boundary z ' = ae as 



_d|_ _ 2d E _1_ 
dz' rr D z' 



(31+) 



v/here 
E = 



cos^P cosh^Yo (^^^^^-l' " cos'^s) + sin'^-p slnh^ Yo('^°^^^^ " sin^o) 



- i sin 2P sinh 2y^ sin 2 6 



1/2 



D = |cosh 2yq - cosh 2(\i/ + i9)| 

= |(cosh 2Yq - cosh 2\i/.cos 29)^ + (sinh 2\l; sin 29 j^]''"^^ 

Finally, combining in equation (21) the factors 
given in equations (2[|.), (32), (33b), and (5l|.) yields 



V 



= ABCD ^ V 
E 



(35) 



v/here 



A = ~- sinh Yi si^ a(cos ^ - cos o^) 

cosh 2y-j^ - co.'3 20 I- 

+ cosh Yi CO.: a (sin CJ - sin 



1? 



C = 



2 2" -1/2 



D"[^co3h - cosh 2^ ' cos 29)^ + (sinh 2^ sin 28)^] 

E= cos^P cosh^YoC^^*^'^'^^ cob^Q) + sln^fi isinh^Y^ (^cosh^xl^ ^ sin^e) 



- ^ sin 2P sinh 2Y^ Qln'2G 



An application of. equation (55) ^'or the purpose qf 
illustrat5-ng the varrious steps involved in a calculation 
of the surface velocity and: pressure of .the airfoil lattice 
is given in appendix E and illustrated in figures 5 6. 
For the sake of coinparlson, the single-airfoil- ■ case is 
given in figure ?• 

Some spec i al res ults froivi equation (35) ^ lattice 

of line s in the case of a lattice of stx'^aight lines, the 
z^-, and w-planes merge* hence 9 = cp =: a and 

R = S = a. " ' - ■ ■ 

Prom equations. ( I9 ) ^^-nd (7) and with a + P '= a^., " 
which is the angle of attack with respect to the chord, 

p sin a.- 

i:~T::ir^ ^^^^ 



2Vd (aosh^Yo - sin^s) 



The lift per unit span on a single raeraber of the lattice 
is given by 

L - pvr 

where p is the air density. The lift vector is per- 
pendicular to the mean velocity vector (fig. 1;.). This 
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result ia general and not limited to a straight-line 
lattice. The lift coefficient is 

pvr 2r , 1 r , 

Cl = --r r = — = I}. — ■ — (57) 

e^|pv2j c/d2Vd 

v/here r/2Vd is given in equation (56) and c/d can be 
found by equation (8), 

The local velocity on the surface (equation (55)) 
becomes 

V = V^os a^^ + I sin a^^. (38) 

where 

jj » ^ cos 6 cos 9 ^ sin P sin 9 

(cosh^Y^ - sin^e)^^^ Yo 3inh Yo 

COS 3 sin 8 sin 3 cos 9 
M = — • 

sinh Yq cosh Yq 

In the special cases in v/hlch P ='0^ and P = 90^, the 
relations (56) to (58) are simpler. 

For stagger angle P = 0*^ and with d = 

p ^ sin 
2Vg cosh Yq- 

From equation (8), 

rrc 



and 



cosh Yr^ = coth 

^ 2g 



L = 2py^g tanh — sin a-_ • 
2g ^ 

tanh 51 

= TrpcV"^ ^ sin 

rrc 
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The lift coefficient, according to equation (5?)* is 

tanh 

Ct = 2tt ~ sin a„ 

TTC 

o 

For p = 0 , therefore, the -slope of the lift curve is- 
alv\fays less tlian 2ir. ..Note that, for large gap,, c/g — ? 0 
and uhe lift coefficient is 

Cl = 2Tr sin. . . 

When the gap g is snail compared v¥lth the chord ^ c, 

' ■ . Ct ]i S sin 
c 

The local velocity at the surface, by equation (38) ^ is 

V = v(cos Oy^- + tanh cot |- sin a^^j 

This result may be conipared with that for the sirlgle-line 
airfoil ( y^. 

, 9 . \ 

+ cot r- sm CL. ) 

2 y 



V = V |cos a.^ 
For stagger angle (3 = ^0^ and with . d- T^ h. 



sm a^- 



^ 2rn sihh y 

: . ■ . VO . 

Prom equation (8), 

• o 2ii 



and 



L = 2pV^h tan ^ sin 



tan HP- 
- Trp®V- ■ ■ sin ax 

2h 
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The lift coefficient, according to equation (57)^ is 

TTC 



tan 



TTC 
2h 



2h : 

sm 0^ 



Por P = J therefore., the slope of the lift curve is 
always greater than 2Tr. The local velocity at the 
surface is 



V = V 



cot 



ay^ + coth Yo cot - 



sm a 



'X 



It may be noted in passing that, for c = —h. 



as compared with 



Cj^ = 8 sin 



= 2tt sin 



for the single airfoil. 

For the limiting case in which b and d approach co;; 
the transformation (6) becomes 



2Trb 



z ' e 



d -13 
and, with limit — i and a new variable z'^=z'e 

2TTb 



a*" 



which is the familiar Joukowski transformation. If the 
variables if and 9 are introduced, the corresponding 
result is expressed as 



i = 2a cosh + 1(6 - B)] 
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where the llirdt, as Yq — ^ "^^ been 

2rra cosh Yq 

put equal to 1- 

to • 
to 

Langley Memorial Aeronautical Laboratory, 

National Advisory Gonurdttee for Aeronautics, 
^Langlejr Field, Va. 
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APPENDIX A 
mil} SYI.©OLS 



i complex plane of airfoil lattice (x + xy) 

i-^, ^2 complex planes of airfoil lattice for stagger 
angles p = 0° and p = 90°, respectively 

(^1 ^yi 5 ^2 ^^2) 

ae / 

z circle of radius R = ae " in z-plane Vae " / 

w circle of radius S = be"^l in w-plane 

^ =co, =b, z= P^, w=:b corresponding points 

i = -00, z ' = -rb, z = -^2' w = -b corresponding points 

a, b reference lengths 

Yq gap-chord parameter (b =: ae^*^} 

P stagger angle 

d lattice spacing, or *'slant" gap for any value of P 

g lattice spacing,, or gap for 3 = 0^ 

h lattice spacing, or stagger for p = 90^ 

V magnitude of mean of inlet- and outlet-velocity 

vectors (fig. ii) 

angle of attack with respect to x-axis of mean 
velocity vector 

a angle of attack with respect to normal to slant 

line of lattice of mean velocity vector 

^1' ^2 inlet and outlet angles of attack v/lth respect 

to normal to slant line of lattice, respectively 

V]_, V2 magnitudes of inlet and outlet velocities, 
respectively 
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APPENDIX B 
INVERSIOII OP EQUATIONS (!•) TO (6) 

AND CHOICE OF COORDILIATES 



It is desired to find from a. given airfoil lattice 
in the ^-plane the contoiir defined by ^1^(9) in the 
z ^ -plane . This problem corresponds to an inversion of 
equations to (6) and can be exactly treated for the 

cases in which ^ - 0^ and P = (equations (li) 

and (5)^ respectively) but an iteration or successive- 
approximation method is required for equation (6). 
Furthermore, although the parametex^s g and h are 
fixed by the geometry of the lattice, a choice exists 
in the definition of the chords and the origin of coordi- 
nates. This choice is discussed follov^ring equation (317) 

Stagger ang le P - 0^ . - Prom equation (5), there is 

obtained 



cosh + l6) = cosh taniri — ^-i 

o 6 



(Bl) 



Putting . t-^ -^1 ^^1 ^^^^^ denoting the real and imaginary 
parts of equation (Bl) by ^-^ and r]-,, respectively, 
leads to 



cosh ^ cos 9 = 



sinh \1/ sin 6 = T)^ = 



2Tr 



cosh Yq sinh — x^ 



g 



cosh — X-, + cos - — 



S 



2ir 



cosh Yo sin ~ y-L 



^ 2Tr ^ 2rr 
cosh — + cos — 

g 1 S ^1 



y (B2) 



The expressions containing x-i^ ar^d y-| in equation (B2) 

are considered given since these quantities are known 
from the coordinates of the airf oil lattice • If \|/ 
and 9 are eliminated successively, 



21+ 



5l V 



^cos Qj 



^1 
,sin 8/ 



= 1 



and 



> 



\cosh W \sinh ^/.■ 



(B5) 



Prom equation (B5), there result the follov/ing expres- 
sions, which serve to define the function ^K9 ) in terms 
of the airfoil coordinates: 



sin^G = p + \/p^ r\^^ 



> 



where 



1/. ? 2 2\ 



(Bl|) 



For small values of 0,., the relation slnh = 
mav be used. 



= ^1 ■ 



sm 



It is useful' for coriiputational purposes to record 
the real and imaginary parts of equation (5) 

2"" 



3. I 1.. £l 



^1 ~ 2TI- \ 2 



log 



> 



^^1 = ^ (^1 - ^2) 



(B5) 



where 



(cosh Yq cosh ^ cos 6 )^ + (sinh sin 9)^ 
(cosh Yq cosh \|/ cos 9)^ + (sinh \|/ sin 9)^ 
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sin cp, = — sinh ^ sin 6 
1 Pi 



sin cpp = - -i- sinh \|/ sin 9 
. P2 



The angles are to be chosen betvveen -tt and tr, and 
the quadrants may be determined by noting also the 
relations 

cos cp^ = (cosh Yo ^osh \l/ cos 9) 



cos Cpp = — (cosh - cosh ^1; cos 9) 

Stagger angle p :r QQ^ ,- From equation (S), there 
is obtained 



linli (xl/ + 19) =z sinh y^ tan U 

o h 



(b6) 



With ^2 " ^2 ^^2 ^^-^ real and imaginary parts of 

equation (B6) denoted by ^2 Ti2^ respectively. 



sinh ^ cos 9 = £^2 = 



cosh 1^ sin 9 = 



T)2 



sinh Yo sin ~ xA 
Q h ^ 

cosh — yo + cos -rr 

h ^ h ^ 

sinh Yo siiih — y2 



T- Sir J 2tt 
cosh 72 cos -j^ ^2 



y (B7) 



If \|f and 8 are elir.iinated successively. 
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( 



M - (^i ■ 

cos 0/ \.sin 6/ 
\sinh ^/ \cosh \|ir/ 



-1 



(b8) 



Prom equations (b8) there result finally the following 
expressions, which serve to define the function ^(6) 
in terms of the airfoil coordinates: 



cos^e = q + ^2 + 



(B9) 



v/here 



-Ik 



n ^ 2 2 

1 - ^2 " '^^a 



For values of 9 near ±90^, the relation sinh ^ = 
may be used. 



cos 9 



It is useful for computational purposes to v/rite the 
real and imaginary parts of equation (5) 



Jz =- 



where 




> 



(BIO) 



V 



(sinh Yq + sinh cos 9)^ + (cosh \|/ sin 9)^ 
(sinh Yo ~ sinh ^ cos 0)^ + (cosh \|/ sin 9)^ 
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sin = — cosh ^ sin 9 
5 P3 



sin CP, 



cosh sin 6 



The angles are to be chosen between * -rr and tt, and the 
quadrants, may be determined by noting also the relations 



cos 



(sirJ:i Yo sinh \i/ cos 9) 



cos cpi = ~« (sinh Yq slnh ^{/ cos 6) 

k 



. Arbitrary stag ger angle P and c hoice of c oordi- 
nates > - Because of the traiiscendental nature of equa- 
tion a direct inversion expression seems unobtainable: 
hov/ever, the values {'^\f , 6) that correspond to coordinates 
(x, y) may be obtained v/ithout difficulty by an iterative 
process. For this purpose and for the purpose of choosing 
the coordinate axes, expansions of Xj^, y^^, and y^" 

in powers of ^ are useful. The follov/ing expansions 
may be readily verified: 



d 'o 
x-i « ~ cos : p 



log 



cosh Yo ^ 
cosh Yq - cos 8 



+ ^^ cosh Yq ^'^^ ® 



sinh'^Yo " sin^^Q 
^cosh^YQ - cos^e} 



(Blla) 



X2 - ^ sin P 



sin 9 

2 tan ^ 

sinh Y>^ 



P 2n 

2 cosh'^Yr. ^ 
+ sinh sin 6 ^ . 

/ P 2.\' 
ycosh Yq ■ ^ ® y 



(Bllb) 
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Then 



^ d « 2 cosh Y_ sin 6 . , , 

y ~ ^ cos P ^2 ^ (Bllc) 

cosh Yo ~ cos^9 



J 2 sinh Y/^ cos 9 
y » - d gin e _Io Blld) 

cosh cos 8 



7 =^ Yl 72 

-|^^/P(e) (B12) 



where 



cosh Yo P ^ " sinh Xo sin P cos 9 

P(6 ) =^ 7 Z — ' 

p o 
cosh Yq - cos'-6 

If the x-coordinate of the straight-line lattice, v/hich 

is considered the skeleton of the airfoil lattice^ is 
denoted by Xq, then- Xq is given by the value of 

X = x-L + X2 for = 0, or 

d / cosh Yn cos 6 .1 sixi 0 \ ■ 

XQ-^ (cos p log +2 sin p tan ^ -iifL^ ) (B15) 

^'^ \ cosh Yo ■ cos 9 sinh y^J 

and 

^ - ^© 2I? GO) ■ (BlL). 

where 

cosh cos 6 cos 9 + sinh sin B sin 9 
G(9 ) = 2 ^ _^ i 

cosh Y • cos 9 
'o 

(cosh Yo cos- p sin 9 • sin Yq ^in f. cos 8) 2 sin 9 cos 9 

^cosh^Yo • 005^9} 

-F»(e) 
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In particular, the leading- and traillng-edge points 
X = Xj and x = x^, are determined by the values of 

9 = 6-, and 9 = 6^ that may be obtained from equa- 
l Lr ■ 

tions (7b)* Then, 

•'^^ ^0, + I? % (^^5) 



where 



sinh cosh 

and Xq denotes the leading edge of the line given by 
\|/ = 0. A similar expression holds for x^ . 
From equat5.on (B12), for constant ^{/, 

= I W(9) ■ . 

In the neighborhood of the leading edge, therefore. 

For Xq near Xg , there is obtained fron equation (B15), 



1 + _ Aq ^ • • • 



where the follov/ins relations are found to hold: 



TT L iQ-Q 



I 



Ol ItL 'Jg_g TT 0 
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Hence, 

""0 « ^O j "1;^. ^oC^ 

Then, from equation (Bli|), 



X - 



Xq « Xq - XQ. 



2rr 



2Tr 



-(9 - 9,)' 



G(9) 



0 



It follows from equation (Bl6) that, for. x = Xg^/ 



and 



d ,2 ^ 



With this value of yq and equation (BI5), 



^0. 



- 2 



^0 



I 



If the total ordinate for both upper and lower sides at 
X = Xq^ is denoted by y^. 



^1 " ^0- 



(BI7) 



This result leads to a siraple 
choosing axes of coordinates 
behave smoothly at the edges; 
of \|/ at the leading edge is 
the values of \}/ at nearby o 
lov;er surfaces. For a parab 
ordinate through the focus, i 
from the vertex to the focus, 
that the end point of the ske 
mately the focus of a parabol 



and convenient v^^ay of 
in order that \l/(0) will 



thai 



that the value 



approximately the mean of 
rdlnates on the upper and 
ola the latus rectum, or 
s four times the distance 

Equation (BI7) states 
leton chord should be approxi- 
a at the nose. 
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The scheme for choice of axes is as follov/s: Locate 
a point P near the leading edge where the ordinate 
through P is four times the distance of F from the 
leading edge. Similarly locate a point P' near the 
trailing edge. The origin of coordinates then bisects 
the line PP^ which is on the x-axis and represents the 
chord of the skeleton line airfoil \|/ = 0 . (To the order 
of approximation employed, the afore-mentioned choice of 
axes coincides with that given for the single v/ing section 
in reference 1 or y.) 

Procedure for finding 9 ) f rom (x, y) for 

arbitrary stagger angle p.- An iterative procedure is 

given herein for finding \j/(9) from (x, y) for arbi- 
trary in w^hich the knowledge of the case for 3=0^ 
is employed to help in formulating the initial approxi- 
mation. - In brief, values of 9 are obtained for ^s tagger 
angle p = 0^ for both the airfoil and. its line skeleton. 
Values of 6 are then found for the skeleton, in the case 
of stagger angle 6. These functions permit approximate 
values of 9 to be found for the airfoil, for stagger 

•angle P# Equation (B12) then enables approximate values 
of ^ to be obtained. These values of (\!/, 6 ) are then 

"readily checked and improved, if necessary. The steps 
are as follov/s : 

(1) Choose the axes as outlined and express the air- 
foil coordinates in percent chord, where the chord for 
this purpose is the part of the x-axis intercepted by the 
airfoil. Denote the coordinates thus obtained by • 

*(xp, yp) . Find, k .= PP» in percent chord. Pind. 

xj - ^07^ distance from the leading edge to P ■ in 

percent chord, and denote. this yalue by e. Obtain the 
ratio c/d, where, c means here PP' and d is the 
spacing between corresponding points on adjacent airfoils 
of the lattice. Find conversion factor m by 

* m = 2ir 5 f . 

d k 

(2) Convert coordinates of the airfoil from /xp, yp^ 
to (2tt |, follows: 
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2lT 2 = my 



(5) Find the parameter Yq that corresponds to 

the determined value of c/d for the given value of P 
from graph or by calculation ( equation (3)), Also find 
for later use the value of c/g corresponding to this 
value of Yo • p = 0^. 

([|.) Consider, for this value of Yq^ ^^^^ 
straigtit-line cases {^\f = 0, 3 = 0^) and = 0, 6 = pj; 
associate values of -9 = 8q for P = 0^ v/ith values 

6 = Bp for the stagger angle P by referring associated 

values of 8 to geometrically similar points of the lines 
(equation (BlJ)). 

(5) Multiply coordinates in step (2) by the ratio 

T where the chord-gap values are from step (3) for 

(c/d)p 

P = 0^ and for P P . Using equation (B^), find 

values of 6 for p = 0°. 

(6) Vvith the aid of step (ij.), obtain approximate 
values of 9^ associated v;ith the values of 9 obtained 

in step (5)- Then, with 8 = Bp, use equation (B12) 

to obtain an approximate value of \|/, vmere 



Zvy F(e) 



d 2 



and the leading- and traillng-edge values of \l/ are 
obtained from equation (B15)* 

(7) Calculate, from equations (B5) and (BIO), exact 
values of — , Zrr ^J, associated vvith the initial 

values of (\!/, 9) iii step (6) where x = x-^ + X2 and 
y = 71 + 72- 
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(8) If^ on comparison- of the coordinates in step (7) 
with the coordinates in step (2), it is deened necessary 
to approximate (.xj/, 0 ) more closely for several of the 
points', (x, y), one procedure is the following: An ex- 

— ^ can be found from equations (Ij.) 



pression for 
to (6) as 



d{^\f + 10 ) 



d(\j/ + i9) 2tt 



sinh (\};+ i9) ^ sinh (\|/ + iS ) 



cosh Yq + cosh (\|/ + i9) cosh Yo ~ cosh + 10) 



cosh {^if + i8) 



d i cosh ( \l/ + 19) 

_ ^ -sin 3'-~ ■ — .- 

2!T jsinh Yo **• sinli + 10) sinh Yq" ^^^^ ii^ + i9) 

With the notation of equations (B5) ^nd (BIO), this ex- 
pression may be vvritten 



d(\i/ + ie) 



= p + iq 



COS 8 sinh + 19 )( -f- e 



-ICP- 



1 -'-"^2 



e 



- 1 sin 3 cosh + 19 ) 



where 

p = cos ^ 



1 .-i*5 . 1 e-n' 

'5 "k 



sinh ^ cos 



+. cosh. ij/ sin 9 



/cos cPt cos cp^\ 
9 ( i + d 

V Pi 



'sin cp^ sin 



+ sin P 



"a mil * s m 6 {22111 * '-^Llk' 



- cosh \k cos 



3k 



and 



q = cos .6 



cosh \|/ sin 8 



/cos cp^ cos cp^X 



4 w . n /sin (P-1 sin cpp'' 
- sinh ^ cos 9 f + 



- sin P 



P2 



/cos Cp7 COS cp].'" 
cosh \1/ cos 9 f - ? + '± 



^ /sin ^>7.■ . sin Cp),N 
+ sinh \j/ sin 9 ( — ■ ^ + ^ 



^5 



The following relation may then be- noted 
+ LA9 » 



p + iq 



Let 



(B18) 




where the subscripts 0 and 1 refer to the coordinates 
given in steps (2) and (?), respectively. If the values 
9) obtained in step (6) are used, evaluation of equa- 
tion (Bl8) gives values AG), and + , 9 + A9) 
represents the next approximation to the desired coordi- 
nates. The process in steps (7) and (8) can be repeated 
if deemed necessary. 



APPENDIX C 
DERIVATION OP EQUATION (12) 



The transformation (equation (9)) from the z«- to 
the 2 -plane, may be rewritten 



log-- 



= f (z) 



(CI) 



where the complex constants c^ may be defined as 



Cv> = a„ + ib^ 
n n n 



On the boundaries, z' = ae 



>ko+lcp 

and z = ae ; hence. 



log 



at 



^ - ■^Q + i( e - cp) 



ana 



_ \ 



/5i 



n 



(C2) 



v/here 



R = ae'-o 



With ^- considered as a fujiction of <p denoted by 
'•i/(cp), the coefficients in equation (02) are obtained as 

;>2Tr 







Rn ~ 




°n _ 


1. 


R" ~ 








R^ 


Tf 



\[r(cp) cos nqp dcp 



.2nr 



n2Tr 



\i/(®) sin ncp d(p 
xj/dp } e a<ip 



(03) 
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Substituting equation (C5) in equation (CI) yields 

f(z) = i > ^ ^„ dq) (Ci^) 

< 1, the geometric series in equation (Cij.) can 



For 



R 
z 

be suhmied and 



^2Tr ^^i<P 



f(z) = i T"" ^{r(cp) . dcp (C5) 

which can immediately be expressed as in equation (12). 

Pop computational purposes, equation (12) may be 
separated into real and imaginary parts. Let 
f(z) = p + Iq and z = x + ly (where, for example. In 
the zeroth approximation x = b, y = 0). Then, 

1 p2Tr K 
P = - *(cp) i<, 



i r 



2tt it, 



p 



\l/(cp) drp 



v/here 



Jo D 



Ni = 5 COS cp + I sin cp - 1 

R R 

No = - sin cp - 2 cos cp 

R R 

D = 1 - 2 Z' ^ cos cp + 2 sin cp'^ + '*'J^ 

\R R / p2 



and the integrations can be conveniently performed by 
Simpson's rule. 
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■ APPENDIX D 
TRANSFORMATION PROM z-PLAI:IE TO w-PLANE 

The linear fracti6nal transformation 

az + b 



V/ 



C2 + d 



on which the derivation of equation (I3) is based, 'has 
the following v/ell-known properties: 

(1) When z traverses a circle Cg, w traverses 
a circle C^^. 

(2) Two points w^ and ^2 inverse with respect 
to a circle correspond to two points z^^ and Z2 
inverse with respect to the circle • 

(5) The anharmonic ratio of four points is pre- 
served; that is, if z-L, Zj, and Z|^ correspond to 

I, Wj-, and wi^. 



w 



(^1 ^ ^h){^2 ^ ^-5) ^ frl ■ ^^^U)fr2 - ^^3) 
(zi - Z2)(25 " ^k) (^'1 " ^'2)(^^3 - 

Ppr the desired correspondence it is knovm that four 
points w^ = b, v/p = -b, and their inverse points 

= — , Wi = are to correspond to Zn = Pi , 

5 b 4 b 11^ 

p2 p2 
Z2 = -P2 ^nd their inverse points z^ = ~, zj = 

^1 ^ ^2 

Property (J) yields a relation that may be used to solve 
for the radius S and that can be expressed by equa- 
tions CI5) and (l6)* Vi/hen the radius of the circle in - 
the w-plane has been determined, property (5) can again 
be used by replacing - say, wj^ hj w and by z. 

This procedure will 'yield a result that is equivalent to 
equation (13 ) • 
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APPENDIX E 
OUTLINE OP CALCULATION PROCEDURE 



(1) List airfoil-section coordinates in percent 
chord. 

(2) Choose axes (appendix B, paragraph follov/ing 
equation (BI7)). 

(3) List stagger angle P and find Yq value 
of c/d for the skeleton line lattice (table I, fig, 2, 
and equation (8) ) . 

(1|) Find (^l/, 6) (appendix B). 

(5) Find €(cp) (equation (10)) by method given in 
appendix of reference 3. 

(6) Plot ^ against cp where cp = 0 + c. Find 
constant (equation (11)) and R = ae^°. 

(7) Find complex constants B.^ and ^2 (equa- 
tion (12) and appendix C), 

. (8) Find constants cosh y-^, Y]_^ and 
K = k^^ + ik2 (equations (I6), (I5), and (lij.)). 

(9) Find coraplex constant 6 =|5!e^'^ (equation (26)) 
and real constant \ (equation (28)), Also obtain 
functions m(cp) and |j,(cp) from equation (29). 

(10) Find a and, in particular, (equation (30)). 

(11) .Evaluate factors B, C, D, and E (equation (55)). 

(12) Evaluate factor A in equation (35), first 
choosing the angle of attack . a • as indicated in the • 
following paragraphs: 

The lift coefficient is as in equation (37) 
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Here c/d refers to the value of x/d at 0-percent chord 
minus x/d at 100-percent chord. By using equation (19) 
for r/2Vd, C-j^ raay be- expressed as 



v/here 



and 



H. = 



r) = 



Cl = H sin( a + T)) 



(El) 



\cosh »f^J ysinh Y]_/ 



1/2 



-1 / sin o,^ cosh ' 

tan--^/ i 

\^co3 sinh 



This relation may be used to find a for any desired 
value of Cj^ and it is further noted that a = -rj is 

the angle of zero lift. 

The *^ideal^^ angle of attack, introduced b7/ Theodorsen, 
is defined for a thin section as the angle of attack for 
which a stagnation point exists not only at the sharp 
trailing edge but also at the sharp leading edge. For 
thick airfoils, the ideal angle of * attack is defined in 
the same manner (the pressure difference at the leading 
edge vanishes) although the point that is considered the 
leading-edge point is not precisely defined. If this 
point is taken to be the intersection of the x-axis with 
the airfoil leading edge, the ideal lift and ideal angle" 
of attack are found as follov;s: Let o be the value 

of o corresponding to the leading-edge point. The 
quantity dQ/dw in equation (or the factor A in 



equation (55)) vanishes for a = a- 



?he relation that 



gives the value' of the ideal angle of attack 
is then 



sin a _ cosh Y]_ sin Q-j^ - sin cr^ 

cos a sinh yt cos - cos 

* 1 1 o 

and the ideal lift coefficient, from equation (El), is 

Cr = -ii ^ i cos ^ (o. ^ o \ 



v/here 

= posh n cos 1(^1 + ^o)l + ^sihh n sin | j^a^ + a^)j 

(15) The surface velocity ratio v/v is now found 
from equation (35) the (mean) superstream pressure is 

found from Bernoulli's equation as 

The angle through v^rhich the strearPi is turned raay be found 
from equation (25). 

A remark may be inserted here regarding an inverse 
calculation procedure. Instead of starting with a given 
lattice, it may be convenient to start with given func- 
tion ^^^(cp), (quantity t as- a function of q)) and given 
parameters and P. Then both the lattice arrange- 

ment and the flow properties follow uniquely, and in this 
way, systematic families of lattices can be studied. 
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